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M. Uplaznik, B. Bercic, M. Remskar and D. Mihailovic 

Jozef Stefan Institute, Jamova 39, SI- 1000, Ljubljana, Slovenia 

Charge transport measurements on flexible M06S3I6 (MoSI) nanowires with different diameters 
in highly imperfect 2-terminal circuits reveal systematic power law behaviour of the conductivity 
cr(T, V) as a function of temperature and voltage. On the basis of measurements on a number of 
circuits we conclude that the behaviour in thin wires can be most convincingly described by tunneling 
through Tomonaga-Luttinger liquid (TLL) segments of MoSI wire, which is in some cases modified 
by environmental Coulomb blockade (ECB). The latter are proposed to arise from deformations or 
imperfections of the MoSI wires, which - in combination with their recognitive terminal sulfur-based 
connectivity properties - might be useful for creating sub-nanometer scale interconnects as well as 
non-linear elements for molecular electronics. 



I. INTRODUCTION 

While the transport properties of one-dimensional sys- 
tems have been of great interest from the point of view 
of fundamental physics for some time^, recently, further 
interest in the transport properties of nano-scale one- 
dimensional systems was aroused because of their im- 
portance for the development of molecular electronics, 
where diverse molecular devices (switches, memory el- 
ements, sensors) all need to be self-assembled together 
with electrically conducting molecular-scale wires. To 
be of practical use, the connectors need to have reliable 
contacts and also be able to withstand mechanical defor- 
mations while retaining their conducting properties. Till 
now there has been no recognized material which could 
be used for this purpose, and this has seriously impeded 
progress in the development of large scale molecular elec- 
tronics in recent years. 

In this paper we investigate the transport properties 
of MogSalg molecular wires^, which have been recently 
shown to be very promising flexible molecular-scale con- 
ductors. MogSalg wires are air-stable one-dimensional 
inorganic cluster polymers (see Fig.l), which are unique 
in that they enable covalent bonding to gold surfaces 
and organic molecules via sulphur atoms at the ends of 
each molecular wireP. Single molecular wires were also 
recently shown to self-assemble gold particles into criti- 
cal scale-free networks^. Such molecular wires may ex- 
pect their electron transport properties to be governed 
by quantum properties on the microscopic level. Thus, 
to make further progress in molecular electronics with 
MoSI connectors, we must first investigate and under- 
stand their molecular scale electronic transport proper- 
ties. Earlier experimental work has shown metallic sig- 
natures, such as a low-frequency Drude response in the 
optical conductivity^, but relatively low room temper- 
ature conductivities a ~ 10 S/m, which decrease with 
decreasing temperature^. Recently Venkataraman, Hong 
and Kim-2 described electron transport measurements in 
multichannel I^MogSeg nanowires with diameters in the 
range 7.2 to 12 nm in terms of a Tomonaga-Luttinger 
liquid (TLL) in contact with Fermi liquid (FL) elec- 
trodes. In another related material NbSe3, nanowires 



ranging between 30 and 300 nm in diameter have been 
observed to display similar power-law behaviour, albeit 
over a rather limited range of temperature^. The extreme 
one-dimensional nature of the MoSI wires suggests that 
signatures of TLL behaviour might also be observed in 
thick multi-strand bundles, and not only in very thin 
wires. On the other hand, we may expect that the elec- 
tronic transport properties in MoSI wires might depart 
from ideal TLL behaviour because of the deformable na- 
ture of the S bridges which link together the M06 clusters 
into ID chains^. 




Figure 1: The Structure of M06S3I6 wires. The axes are in 
the frame of the monoclinic space group P;P^ with center- 
to-center distance of do = 0.958, and the unit cell length 
c = 1.197 nm. The S atoms linking the Mo clusters form 
accordion-like deformable bridges and are shown in yellow, 
while the orange and violet spheres represent I and Mo re- 
spectively. 

In this paper we present a systematic study of MoSI 
wires of different diameter from 4 nm to 1000 nm, ex- 
amining the T-dependence, diameter dependence and 
current- voltage characteristics at different temperatures. 



We are particularly interested in the behaviour of realis- 
tic circuit configurations, with irregular wire geometries. 
The intrinsic flexibility of t he Mo SI wires, arising from 
their accordion-like struct ur d 10 ! 11 ! means that they bend 
easily to conform to surface contours. We have therefore 
focused on dielectrophoretically deposited thin wires over 
contacts in which the wires conform to the surface relief. 

Considering the possible transport mechanisms, we 
confine ourselves to the common ones discussed in 
literature^!, namely TLL tunneling^, environmental 
Coulomb blockade^ (ECB) and variable range hopping 
(VRH) in the presence of Coulomb charging effectsSS 

1. The T— and V— dependence for tunneling into a 
ID TL liquid via Fermi-liquid metal contacts is given by: 
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where a = (g' 1 - l)/4,/3 = (g + g- 1 - 2)/8 and the 
Luttinger parameter g — vp/v p . 7 is a fitting para meter 
that accounts for the voltage drop over the circuitPE^. 
A collapsed diagram of the underlying transport charac- 
teristic is obtained by plotting I/T a+1 against eV/kT, 
where a is the slope of zero voltage conductivity against 
temperature a — a^T a . (3 is the exponent for the high 
voltage limit (eV kT) arising from the power law be- 
havior / oc V /3+1 . 

2. Unfortunately ECB models cannot be solved ana- 
lytically for the general case, but the asymptotic behav- 
ior is very characteristic. The experimentally obtainable 
low-temperature behaviour is given by: 
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where Ec is the charging energy, g — Gq/G, G is the 
frequency-independent conductance, Go = 2e 2 /h and 
Rt is the tunneling resistance.^ For low voltages and 
temperatures the current follows a power law behavior 
/ oc V 2 / 9 . For high voltages (but low temperatures), 
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which gives a linear I — V dependence at high V, so the 
derivative dl/dV is ex pecte d to approach an asymptotic 
constant value of also gives a non-zero inter- 

cept for I — given by the charging energy e 2 /2C. 

3. For the variable range hopping mode, a plot of ln(G) 
(for low voltage) against T~ x yields curves which become 
linear with the correct hopping exponent A. The fits to 
the data typically give a large error in the exponents, 
so to extract the best value of A, fits to the data are 
tested statistically by calculating Pearson's correlation 
coefficients. 

Rather than choosing a few measurements which obey 
one or another type of behaviour, we present here a sum- 
mary of a number of experiments, to try and under- 
stand the different types of behaviour that can arise in 
nanoscale circuits with MoSI wires. 




Figure 2: A schematic diagram of the circuit geometry (left). 
AFM topographic images of the thin wires after annealing. 
The wire diameters (from the height profile) are 4.2 and 4 nm 
for samples na23 (middle) and na27 (right) respectively. The 
roughness of the Ni electrodes appears after annealing and 
is caused by metal aggregation during the high temperature 
annealing process. 



II. EXPERIMENTAL DETAILS 

The thin wires were prepared according to the method 
reported by Nicolosi et al (author?)^ by repeated dis- 
persion and dilution. The dispersion procedure separates 
the wires into two distinct categories: thin wires, with di- 
ameters 4 < D < 10 nm and thick multichannel bundles 
with 100 < D < 1000 nm(author?)lSl. In the ultrasonic 
bath processing procedure, the defective wires break up 
into shorter segments, leaving less defective long thin 
wires in solution. The individual strands within the thin 
wires may thus be expected to have significantly fewer im- 
perfections than within the thick bundles, which may re- 
sult in different electron transport behaviour in thick and 
thin wires. This should be evident in the room temper- 
ature conductivity (<73ooif) as well as the T-dependence 
and I — V systematics. 

The wires were dielectrically deposited typically over 
nickel electrodes prepared by electron beam lithography 
(EBL) by placing a drop of solution over the electrodes 
and applying a 50Hz AC electric field to the electrodes. 
The entire circuits were then annealed in vacuum at 700 
C for an hour. The circuits after annealing are shown in 
Figure 2. Typical resistances of the nanowires at room 
temperature were between 100kf2 and lOOMfL Care was 
taken to ensure good thermal contact of the sample with 
the cryostat cold finger. 



III. EXPERIMENTAL RESULTS AND 
ANALYSIS 

Thin wires 

Two thin wire circuits labeled na27 and na23 (shown in 
Fig. 2) illustrate the most common type of behaviour ob- 
served in a number of measured circuits with diameters 
ranging from d ~ 4 to 15 nm. The room-temperature 
conductivity (T300K for na27 and na23 was 3710 S/m 
and 11900 S/m respectively. Their I — V characteris- 
tics at temperatures between 18 K and 300K are shown 
in Figs. 3a) and 4a) and exhibit qualitatively different 
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Figure 3: a) Normalised I— V characteristics of sample circuit 
labeled na27 (d = 4nm) at different temperatures, where Iq = 
230 nA. b) The conductance as a function of temperature on 
a log-log plot, c) A plot of I/T a+1 vs. eV/kT shows the data 
collapse onto a single curve. 



behaviour. In the case of sample na27, the I-V curves 
show characteristic inverted S-shaped curves whose cur- 
vature is strongly T-dependent. In contrast, circuit na23 
shows a clear J-like shape characteristic. (Other cir- 
cuits we have measured show behaviour in between these 
two extremes^) . The T-dependence of the conductance 
G = dl/dV at V = is shown in Figs. 3b) and 4b). The 
line is a fit to a power law G = GoT a where a = 2.3 and 
3.5 respectively. The data follow the power law fit rea- 
sonably well, but do not give a perfect fit over the entire 
range of T. Plotting the entire data set l/T a+1 against 
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Figure 4: a) Normalised I — V characteristics of sample circuit 
labeled na23 (d — 4.2nm) at different temperatures, where 
Io — 600 nA. The arrow shows the extrapolated voltage cor- 
responding to V — e/2C. b) The conductance as a function of 
temperature on a log-log plot, c) A plot of I /T a+l vs. eV/kT 
do not show the data collapse exhibited by circuit na27 in Fig. 
3c). 



eV/kT according to the TLL prediction, for na27 we see 
that the data collapse quite well onto a single curve (Fig. 
3c), where (3 — 1.6 is obtained from the fit. The be- 
haviour of na23 is very different to na27. It gives no 
such TLL collapse (Fig. 4c), indicating clear departure 
from TLL predictions. (Overall, approximately half our 
circuits showed the TLL collapse.) 

Attempting to fit the ECB model to the data, we 
would expect the I/V slope to cross over from ~ 2/g 
at low voltages to 1/Rt at high voltages. Correspond- 
ingly, the derivative dl/dV should show a systematic T- 
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Figure 5: The derivatives dl/dV obtained from the I-V data 
plotted in Figs. 3a) and 4a) for circuits a) na23 and b) na27 
at different temperatures. The derivative at low V increases 
in both cases, as T increases from 18K to 300K 
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Figure 6: I-V characteristic of sample circuits a) na23 and b) 
na27 at different temperatures. 



independent cross-over from 2/g — 1 to corresponding 
to the low and high voltage limits respectively. We plot 
dl/dV vs. V, in Figs. 5 a) and b) respectively for na23 
and na27. Circuit na27 shows a rather small variation of 
dl/dV with V over the entire range of V and T (note the 
small values on the log scale on the dl/dV axis), quite un- 
like the ECB model predictions. In contrast, na23 shows 
a systematic variation with V, which for low tempera- 
tures is not far from what is expected on the basis of 
ECB theory: at low V, dl/dV ~ (within experimental 
error) , but clearly saturating at high V. 

For completeness, attempting to fit the data to the 
VRH model, we first plotted ln(G) (for low voltage) 
against T~ A for different A=l/4 to 1, according to the 
VRH model (Fig.6). Making a statistical evaluation us- 
ing Pearson's correlation test, we find that the best fit is 
with A = 1/4 for na27 and A = 1/3 for na23. Such small 
values of A within VRH models would imply the dom- 
inance of 3D hopping processes, which is not expected. 
(We shall discuss the implications later). 



Thick bundles 



A ubiquitous feature of the thick wire bundles (d = 
100 ~ 1000) is their linear I-V characteristic from 18K 
to 300K®. The room temperature conductivity <jq is typ- 
ically 3 orders of magnitude smaller than for the thin 
wires, around 1-10 S/m. This is taken as a clear indica- 
tion that only a small fraction / of the molecular strands 
in the bundle contribute to the transport. Confirming 
earlier preliminary measurements, the T-dependence was 
not very well described by a ID VRH modeP. The new 
systematic data on many circuits now confirms this. The 
dependence of the conductivity is shown in Figure 7 for 
a number of bundles of different diameter. Surprisingly, 
the data for the different bundles all appear to follow 
power-law behaviour cr = aoT a quite well. Moreover, a 
systematic trend is observed, whereby the larger diam- 
eters have a smaller exponent a, an indication that the 
numbers of conducting channels scales with diameter^ 
(Figure 8). 
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Figure 7: The normalized conductance G/Go of a large num- 
ber of thick wires, with d = 200 ~ 2000 nm. In two mea- 
surements, a break is obseved, which is believed to be strain- 
induced change in the number of conducting channels. 
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Figure 8: The dependence of the temperature exponent a of 
thick wires on bundle diameter d from the data in Fig. 7. 



IV. DISCUSSION 

The present set of experiments shows a rather wide 
range of behaviour, suggesting caution in trying to inter- 
pret the data by any single model. Discussing the simple 
VRH model first, our fits give 3D exponents, which is 
clearly inconsistent with the nature of the system investi- 
gated here and is dismissed from further discussion. The 
calculation of Fogler et al.^ considered a quasi-lD wire 
with a finite density of impurities modeled by a series of 
weakly coupled quantum dots. They predict power law 
behaviour of the current for both T and V, as / cx V@ +1 
at high V and I cx T a V at low V with a > /3 > 1. The 



model holds for a large number of parallel statistically 
independent conduction channels, where power laws are 
not obscured by the fluctuations of G. This case is ful- 
filled here for thick wires, with bundle diameters between 
100 nm and 1 /i might have 10 4 to 10 6 conducting chan- 
nels. However, the prediction that (3 ^S> 1 is not fulfilled 
here for either the thin or the thick multichannel bun- 
dles, suggesting that this model may not be applicable 
for describing MoSI circuits. 

In circuit na 27 the TLL prediction can be confirmed 
with some degree of confidence by the collapsed plot of 
j jrpa+i vg _ e y/kT. For tunneling from a Fermi liquid 
into a perfect TLL without defects, the model predic- 
tion is that a — (3. Imperfections, such as deformations 
and kinks or stoichiometric defects break up the wire 
into TLL segments, which introduces TLL-TLL tunnel- 
ing between these segments^, for which the predicted ra- 
tio changes to a = 2(3. The measured value for na27 is 
a/(3 — 1.45 ± 0.1. A similar circuit with a wire diameter 
of 5 nm (nal2, not shown) gives a near-identical collapse 
as na27 with a/ (3 = 2.00 ± 0.1. 

The exponent a for a TV-channel wire is given by TLL 
theory®; as: 

where g is the electron-electron interaction parameter. 
The number of conducting channels N in a nanowire of 
diameter d can be written as N — 2/^-, where a is the 
lattice constant, and / is the fraction of molecular wires 
which actually carry current without interruption. For 
circuit na27, using / = 1, a = 2.3, and N = 16 channels 
(for d = 4 nm), we calculate g = 0.21. This is slightly 
larger than the value 0.15 obtained by Venkataraman et 
al for Li2Mo6Se6 nanowires. Other wires which we have 
measured have a between 2 and 3.5, giving a range of 
0.09 < g < 0.27. Comparing this with g estimated from 
materials parameters, we can use the expression from the 
Coulomb charge screening modepS g 2 = vfRqCq, where 
vp = 7 x 10 5 m/s is the Fermi velocity calculated by DFT 
calculations, Rq is the resistance quantum and the single 
wire capacitance per unit length is Cq — ifpf/fb" - Using 
e r = 3.9 for the Si oxide layer of thickness t = 600nm, 
we obtain g — 0.49. In this estimate, the DFT value 
of vp is likely to be overestimated, which may account 
for a large part of the discrepancy between the values 
obtained directly from the measurements. It has also 
been noted previously that the static Coulomb screening 
model overestimates g^. 

The data in Figure 8 for the thick wires can be fit using 
the expression (4). A range of values of g can be obtained 
from the fits to the present data, 0.1 < g < 0.5 with 
corresponding filling factors 10~ 5 < f < 10~ 4 . Values 
outside this range of g cannot be made to fit the data. 
However, since the I — V characteristics are ubiquitously 
linear, for the thick wires, this in itself cannot be taken 
as proof for TLL behaviour. 
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The low value of Co in thick wires is assumed to be re- 
lated to imperfections on the wires. It suggests that the 
transport along the wires is dominated by tunneling be- 
tween TLL segments. In fact the resistance of both thick 
and thin wires is typically of the same order of magni- 
tude, which suggests that in the thick bundles only a 
very small number of wires are uninterrupted, or that 
only strands on the outside of the bundle are conduct- 
ing, with the inner strands being unreachable due to the 
extremely small perpendicular inter-molecular hopping 
rates within each wire. 

Turning to the case of na23, which shows clear char- 
acteristics of ECB behaviour, the intercept of the I — V 
curve at / = gives V — ^ — 0.15V and consequently 
C ~ 5.3 x 1CT 19 F. Comparing this QDOT capacitance 
with an estimate of C for a single wire over a Si ground 
plane with a 600 nm SiC>2 insulating layer and length 
L = 265 nm between contacts, we have C ~ 8 x 10~ 18 , 
for na23. Depending on what we assume for the QDOT 
shape, its size appears to be some fraction (~ 1/10) of 
the distance between the electrodes. 

Considering the possible origin of the ECB behaviour 
and the observed departures from the predicted TLL 
characteristics, we can imagine imperfections and breaks 
of continuity in the wires of diverse origin. The most ob- 
vious and unavoidable effect arises from the deformation 
of the nanowires adapting to the relief of the contacts 
(Figure 2). Their inherent flexibility allows for a sub- 
stantial deformation, which is accompanied by signifi cant 
changes in electronic structure near the Fermi energy^EiD. 
We can thus envisage that the deformations can lead to 
breaks in the continuity of individual channels, and as a 
result the formation of a QDOT in between the breaks. 
In this case the ECB capacitance would scale with the 
distance between electrodes. Another effect, which might 
be important arises from the intrinsic tendency of the 
wires to form discontinuities in the structure, as shown 
in the HRTEM image in Figure 5 a). Clear stripes are 
sometimes observed across the wire bundles, which arise 
from from stacking faults which appear during the growth 
process. Indeed such compositional ordering has been 
recently theoretically predicted^. HRTEM diffraction 



analysis shows that the structure of the wire is identical 
on both sides of the fault, but clearly continuity is bro- 
ken at these points. Sections in between the faults may 
thus act as QDOTs, leading to the ECB behaviour we ob- 
serve. Considering that the distance between faults may 
be a few tens of nanometers, the comparison of QDOT 
capacitance with the static capacitance calculated in the 
previous paragraph suggests this effect may also be im- 
portant. 



V. CONCLUSION 

One of our objectives has been to determine how im- 
perfect, bent and deformed MoSI wires might behave 
in molecular-scale circuits such as may form upon self- 
assembly^. It is clear from the present experiments that 
quantum transport dominates their behaviour. The TLL 
model with its characteristic data collapse of the I — V 
characteristics at different temperatures appears to hold 
well for a significant proportion of the thin wire circuits. 
At the same time clear signatures of the ECB predicted 
J-shaped I — V curves are also occasionally observed. 
Discontinuities in the wires either as a result of bending, 
and/or structural stacking faults within the wire are be- 
lieved to cause the formation of QDOTs, which lead to 
the occasional occurrence of ECB behaviour. The diam- 
eter of the present MoSI wire bundles is small enough 
to allow covalent S bonding to individual molecules^, so 
for the construction of molecular scale circuits, where 
thin and relatively short wires are of interest, they may 
potentially revolutionize molecular electronics. A point 
of interest is the possibility of making variable sizes of 
QDOTs with MoSI wires by stretching them over ap- 
propriately sized regular topological features to produce 
arrays of QDOTs. The other possibility of creating single 
QDOTs by using the tip of an atomic force microscope, 
was already recently demonstrated^. 

We wish to thank D.Vengust for providing samples 
of MoSI nanowires and J.Strle for proof reading the 
manuscript. 
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